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ABSTRACT 

We discuss a possible scenario for radiation mechanism of pulsed GeV 7-rays 
from magnetars. The magnetars have shown frequent X-ray bursts, which 
would be triggered by crust fractures and could release the energy of order of 
Etot ^ 10^^~'*^erg. If the location of the crust cracking of the magnetic field is 
close to the magnetic pole, part of the released energy may excite the Alfven 
wave that can propagate into outer magnetosphere. The oscillation of the 
magnetic field induces the available potential drop S^p ~ 10^'"' Volts, which 
can accelerate the electrons and/or positrons to the Lorentz factor F ~ 10'' 
in the outer magnetosphere. The curvature radiation process at outer magne- 
tosphere can produce GeV 7-rays. If the radiation process is occurred above 
r ~ 5 X 10^ cm from the stellar surface, the emitted GeV 7-rays can escape 
from the pair-creation process with the X-rays and/or the magnetic field. The 
expected luminosity of the GeV emissions is order of ^ 10"^^ erg s~^, and 
the radiation process will last for a temporal scale of years. The expected 
pulse profiles have a broad shape with sometimes sharp peaks. We apply the 
model to AXP IE 2259+586. 

Key words: 



1 INTRODUCTION 

Soft Gamma-ray Repeaters and Anomalous X-ray Pulsars (AXPs) have been discussed in 
terms of "magnetar" model, in which the activities of the neutron star are powered by the 
dissipation of the extremely strong magnetic field (Thompson & Duncan 1995; Thompson, 
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Lyutikov & Kulkarni 2002; Woods & Thompson 2006; Mereghetti 2008). The magnetar's 
emissions mainly appear in the X-ray bands, which are described by the blackbody compo- 
nent (with hard tail) below 10 keV plus a hard power law component above 10 keV (Kuiper 
et al. 2006). The power law components are often explained by the resonant Compton scat- 
tering process of the mildly relativistic electrons and/or positrons below 10 keV (Fernandez 
& Thompson 2007; Rea et al. 2008) and of the relativistic pairs above 10 keV (Beloborodov 
& Thompson 2007; Baring & Harding 2007; Beloborodov 2012), respectively. 

Although no pulsed GeV emissions from the magnetars have been confirmed, it has not 
been conclusive that the magnetars are intrinsically dark in the GeV radiation or present 
sensitivity of the Fermi telescope is not enough to detect any pulsed GeV emission from 
magnetars. It has been argued that the rotation powered activities of the magnetars can 
produce the 7-rays from e.g. the outer gap accelerator (Cheng & Zhang, 2001). However, 
because temperature of the surface X-ray emissions of the magnetars is kT ~ O.SkeV, which 
is much higher than the typical surface temperature kT ~ 0.1 keV of the young pulsars, the 
size of the outer gap, and resultant power of the 7-ray emissions will be relatively smaller 
than those of the canonical 7-ray pulsars, indicating less possibility for the detection of the 
pulsed 7-ray emissions from the outer gap of the magnetar. 

Recently, GeV 7-rays from the supernova remnant (SNR) CTB 109, which will be asso- 
ciated with AXP IE 2259+586, have been found in data of Fermi 7-ray telescope. (Castro 
et al. 2012). Although the emissions from SNR CTB 109 has been suggested (Castro et al. 
2012), the origin has not been confirmed yet. Because the future Fermi observations will 
allow us a more deep search of the GeV emissions from the magnetars, it will be worth 
to discuss the possible mechanism of GeV 7-ray emissions in the magnetosphere of the 
magnetar. 

The magnetar's X-ray/7-ray radiations powered by the magnetic energy have been dis- 
cussed with the emission process near the stellar surface r < 10^ cm (e.g. Beloborodov 
2012). In such a case, the GeV 7-rays cannot escape from the pair-creation processes with 
the X-rays and/or the magnetic field, as we will discuss in section [2l In this paper, therefore, 
we will discuss a possible scenario for the GeV emissions process powered by the magnetic 
energy in the outer magnetosphere. In section O we will discuss the pair-creation processes 
of the GeV 7-rays and critical radial distance above which the produced GeV 7-rays can 
escape from the pair-creation processes. In section [3], we will argue that Alfven wave will 
be excited by crust cracking of the magnetic field and carries the released energy into outer 
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magnetosphere along the background magnetic field line. We will discuss the possible accel- 
eration mechanism due to the propagation of the Alfven wave in the outer magnetosphere. 
We will also estimate the Lorentz factor of the accelerated particles and the typical energy of 
the curvature radiation process. We will calculate the expected pulse profiles in section 13.41 
and 7-ray spectra in section 13.51 In section |U after brief summary of the results, we will 
discuss the possibility of the detection for the pulsed GeV radiation from the magnetars. 



2 PAIR-CREATION PROCESSES AND GEV 7-RAY EMISSION REGION 

The ultra-strong surface magnetic field of the magnetars permits one-photon pair-creation 
process, 7 + i? — t- + e~ (Harding & Lai 2006). The condition of the magnetic pair- 
creation process occurred at a radial distance r from the neutron star may be written down 
as (Ruderman & Sutherland 1975) 

^d(r)singkB ^ nf\/n^m n\ 

^ > X for B,ir)/B, < 0.1, (1) 

where x ~ 0.1, E^^ is photon's energy, -Bd(r) is the global dipole magnetic field at r, and 
5c = 4.4 X G is the critical magnetic field strength. In addition, meC^ is the electron rest 
mass energy and ^kB is the angle between the magnetic field and the propagation direction 
of the photon. Using the dipole magnetic field, we estimate the radial distance, above which 
GeV 7-rays can escape from the magnetic pair-creation process, 

.„^3xl0'(^)-"'(3|^)'7|)"\i,>'/»».ecm.. (2) 

The high-energy photons may be emitted along the magnetic field line by the relativistic 
particles, that is O^b ~ 0, and hence they can initially propagate into the magnetosphere. 
However, because the field lines have a curvature and because the field lines are co-rotating 
with the magnetar, the collision angle increases as the photons propagate from the emission 
point. Under the approximation of the concentric circles of the magnetic field line, the 



collision angle 6kB will develop as sin^^^ ~ s/y i?^ _|_ where s is the propagation distance 
from the emission point and is the curvature radius of the magnetic field line. Hence, 
the photons emitted blow r r^ will be converted into pairs after they propagate s ~ Rc, 
where sin^fc^ cannot be small. 

The GeV 7-rays are also subject to the two photon pair-creation process with background 
soft X-ray field. Using typical observed luminosity Lx ~ 10"^^ erg s^^ and temperature 
kT ~ 0.5 keV (e.g. Kuiper et al. 2006; den Hartog et al. 2008; Enoto et al. 2010), the number 
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density can be estimated as nx ~ Lx/(47rr^c_Ex) ~ 3 x 10^°(Lx/10^^erg s"^)(i?x/0.5keV)"^ 
(r/10^cm)"^cm~^. The optical depth r ~ nxo^^r^ where a^^ ~ (7^/3 is the pair-creation 
cross section with ot being the Thomson cross section, is below unity if GeV 7-rays are 
emitted above the radiation distance 



Equations ([2]) and ([3]) imply that the critical radius above which the produced GeV 7-rays 
can escape from the pair-creation processes is order of Tc ~ 5 x 10'' cm. 

3 EMISSION MODEL 

3.1 The Alfven wave excited by the crust fractures 

Woods et al. (2005) proposed that there are two distinct types of magnetar's bursts, which 
were named Type A and Type B. Type A bursts are frequently seen in SGR bursts, in that 
they are uncorrelated with pulse phase and the energy emitted during primary burst peak 
is larger than the tail energy. Type B bursts, on the other hands, are correlated with the 
pulse phase and the energy of primary burst peak is smaller than the tail energy (see also 
Kaspi 2007; Scholz & Kaspi 2011; Dib et al. 2012). Woods et al. (2005) also speculated that 
the Type A and Type B bursts are triggered by the magnetospheric reconnection (Lyutikov 
2003) and by the crust fractures (Thompson & Duncan 1995), respectively. Furthermore, 
several magnetars sometimes show timing glitches that accompany the radiative outbursts 
(Woods et al. 2004;igdem et al. 2012; Dib et al. 2012; Pons & Rea 2012). In 2002 giant 
glitch of AXP IE 2259+586, for example, the X-ray outburst were consisted of a rapidly 
decay emission in the first few hours with a released energy ~ lO^^ergs and a slow decay 
emission lasting several years with a released energy > 10^^ ergs (Woods et al. 2004; Zhu et 
al. 2008), which are Type B burst. This giant glitch accompanying the radiative outburst 
is likely triggered by the crust fractures, which would simultaneously affect both super-fluid 
core and the magnetosphere, and will inject a energy ~ 10^^~^^erg into the magnetosphere 
(Woods et al. 2004; Pons & Rea 2012). 

In this paper, we consider the energy release caused by the crust cracking of the strong 
magnetic field, i.e. the Type B bursts. We assume that the magnetic field in the crust is 
deformed away from the equilibrium state by an amount of bB. The crust fracture will be 
occurred if the Maxwell stress 5i?-Bcurst/47r, where -Bcrust is the equilibrium magnetic field 




(3) 
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in the crust, exceeds crustal elastic stress. The available magnetic energy density when the 
crust fracture is occurred may be expressed by (Thompson & Duncan 1995) 

where fi is the shear modulus and ey is the yield strain at which the crust cracks. The released 
magnetic energy is therefore estimated as 

At \ / ey Y / Be 



^•"-^'^^"'Hios^j l^j li^j [w^) w 

where i is the typical size of the cracked platelet. Note that with a typical luminosity of 
the observed X-ray emissions Lx ~ 10'^^~^^erg s~^, the emission will last a temporal scale 
of years, r = Efot/Lx ~ 10''~^s . 

The process of the 7-ray emission induced by the neutron star quake was discussed by 
Blaes et al. (1989), who investigated the possibility for the origin of 7-ray bursts. In their 
model, the neutron star quake will excite the oscillation of the magnetic field frozen into 
the star's crust, and the Alfven waves carrying the released energy into the magnetosphere. 
The dissipation of wave energy will accelerate the electrons and/or positrons, which in turn 
radiate 7-rays. 

In the magnetar model, it can be thought that the Alfven wave is excited on the magnetic 
field lines anchored on the cracked platelet. If location of the cracking is close to the magnetic 
pole, a part of the released energy i?tot can be carried by the Alfven waves that propagate 
into the outer magnetosphere r ^ Tc ~ 5 x lO'' cm, where the 7-rays can escape from the 
pair-creation processes. Because the Alfven wave propagates along the closed magnetic field 
lines, the wave may bounce many times between the footprints of the field lines. The time- 
averaged amplitude of the magnetic field corresponding to the Alfven wave energy density, 
~ Etot/Rc may be described as 

/ p \ / R \ -3/2 
SB{R.) 5 X 10^^ G (5) 



1042erg/ V2-10W. 
near the stellar surface. 

As the Alfven wave propagates from the stellar surface into outer magnetosphere, the 
amplitude will evolve as 5B{r) oc A~^^'^{r), where A{r) is the cross section of the oscillating 
magnetic flux tube. Because the background dipole field is proportional to Bd{r) oc A~^{r), 
the ratio of the perturbed and background magnetic fields evolves with the radial distance 
as 
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B,{r) B,{R,) [ B,{r) ) ^ [ B,{r) ) ' ^ > 

In the magnetosphere of magnetar, the speed of the Alfven wave exceeds the speed of 
hght, that is, 

where the p is the mass density. In the limit of va ^ c, we can see that the required current 
density [ia) to support the propagation of the Alfve wave at distance r is in order of (see 
Appendix) 

ia{r) uPUr)mr) _ ^^^^^ / ^ ) ( ^ ( ^^^^^^ ^'^^^A f BdjRsW ^' ^ (g) 



KiGj{r) KVLBd{r) VlO^Hz/ VlHz/ ^ IQ-^ J \ Bd{r 

where icj = ^IB/2'k is the Goldreich- Julian current density, and u is the typical shear 
frequency of the wave, which may be estimated as w ~ Pa-ust/^"^^^^"^ ~ 10^ Hz with pcrust ~ 
5 X lO^^g cm""^ being the typical mass density in the crust. 

We can see that in the limit of va ^ c, the induced electric field is same order of 
magnitude as the perturbed magnetic field, 6E ~ 6B (see Appendix). Hence, the maximum 
magnitude of the induced electric potential due to perturbation of the magnetic field lines 
may be estimated as 

^Si.SB= 1.5 X lO-'^ (m^MM) f (^] Volts, (9) 

where 6i is the displacement of footprints of the oscillating magnetic lines and was estimated 

as 

Bd \ 10-3 J \10^ cm) ^ ^ 

The maximum radiation power may be estimated as 

.,..*......io»(A)(_^)V_l^Vf^)e..-, ,n) 



where I is the total current and Iqj ~ icj 



92 



3.2 Acceleration and GeV emission process 

In Figure [H we illustrate the schematic view of our model for the particle acceleration and 
GeV emissions in the magnetosphere of the magnetar. If the magnetic field cracks the platelet 
near the magnetic axis, the oscillating magnetic flux tube will extend beyond Tc ~ 5 x 10^ cm. 
We will argue that the particles are accelerated by the Alfven wave, which is excited by the 
crust cracking of the magnetic field. It will be possible that the energy dissipation of the 
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Alfven wave accelerates the electrons and/or positrons to a Lorentz factor above F ~ 10''. 
The curvature radiation process of the accelerated particle can produce GeV 7-rays. The 
GeV 7-rays emitted above Tc ~ 5 x 10^ cm will escape from the pair-creation process and 
may be observed as the pulsed emissions. In this section, we will discuss possible particle 
acceleration processes related to the propagation of the Alfven wave. 

First, if the co-rotation motion is ignored Vco = 0, the induced electric field 6E and 
the background magnetic field are a mutually orthogonal and hence no electric field along 
the background magnetic field is induced. However, under the presence of the co-rotation 
motion of the plasmas, such orthogonality breaks down, and the electric field and the current 
along the background magnetic field line are induced (Kojima & Okita, 2004). The typical 
strength of the electric field along the background magnetic field line, 6E^^, is described as 
(see Appendix) 

SEiiir) ^ l3,oSB{r). (12) 

If the Lorentz factor of particle is determined by the balancing the between the electric force 
and the curvature radiation reaction force, we obtain 

^^'^ - [-^^^W) -3.5x10 (^^-^j (^^j [-^J ,(13) 

where Rc is the curvature radius of the dipole field. We can see that the typical energy of 
the curvature photons becomes several GeV, that is, 

47r V3-10V VlO^cmy ^ ^ 

Second, the deficiency of the current (or charge) to support the frozen in condition, 
E ~\- 13 X B = 0, could induce the electric field along the magnetic field (see also Blaes et al. 
1989; Fatuzzo & Melia 1993; Thompson 2006). We can see in equation ([8]) that the required 
current to support the propagation of the Alfven wave easily exceeds the Goldreich- Julian 
current. Furthermore, the typical charge density of the unperturbed particles in the closed 
field region will be coincide with the Goldreich- Julian charge density of the dipole field, that 
is Pe = ■ Bfi/lTic. If the Alfven wave propagates into outer magnetosphere, however, the 
perturbation of the magnetic field line changes the local Goldreich- Julian charge density as 
n-{Ba + 5B) 

pGJ,i = ^ = Pgj,o + ^Pgj (15) 

where Spcj = ■ SB/2n. The electric field along the magnetic field line then will be 
induced, if the primary particles could not supply the required charge density pcj.i- The 
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particles could be accelerated at the "charge-starved" region up to the Lorentz factor F ~ 
10^, such as the particle acceleration process of the 7-ray pulsars (e.g. Wang, Takata & 
Cheng, 2010). The accelerated particles will emit the 7-rays via the curvature radiation 
process. Most of the 7-rays produced below r ~ will be converted into electron and 
positron pairs through the pair-creation processes. Because the GeV 7-rays are emitted along 
the magnetic field lines, the pairs are produced at the convex side of the magnetic field line 
on which primary 7-rays were emitted (c.f. Figured]). The multiplicity is roughly estimated 
as K ~ A X e5^/VmeC^ = 2 x /W^Yo\t){V ■ W)-^ with A ~ rj{b x W)cm ~ 0.1 

being the fractional power spent below r ~ rc. A fraction of the created pairs is probably 
charge-separated, and will produce the required current (|8]) to support the propagation of 
the Alfven wave. 

Finally if the Alfven wave does not efficiently dissipate along the background magnetic 
field lines, the ratio of the amplitude of the Alfven wave to the dipole field will evolve 
as 5B{r) / Bd{r) oc lQ~^{Bd{Rs) / Bd{r)Y/'^ , as equation ([6]) shows. Hence, we find that the 
fractional perturbation of the magnetic field becomes order of unity at the radial distance 
r ~ 10^[5i?(i?s))/10~^i?rf(-Rs)]~^/^ cm. Because the induced electric field is same order of 
magnitude as the perturbed magnetic field, 5E ~ 5B, the total electric field is of order 
of the magnetic field \E\ ~ \B\, where conversion from the electromagnetic energy into 
the particles energy could be possible (Beskin & Rafikov 2000). When the nonlinear term 
becomes to be important, by whatever process, a substantial part of the wave energy is 
probably converted into the electron/positron energy, which in turn radiate the GeV 7-rays. 

In any acceleration processes process discussed above, the typical Lorentz factor will be 
characterized by 

=3x10 (,158^; \-^J \W^) I^To^j 

where we used typical electric field E\\ = 6^p/L with L being the arc length of the acceler- 
ation region along the magnetic field line. 

3.3 Incident X-ray/soft 7-ray emission 

In the present scenario, the relativistic electrons/positrons, which were accelerated in outer 
magnetosphere, will be return to the stellar surface, because they migrate along the closed 
magnetic field lines. Below ~ 5 x 10^ cm, most of 7-rays emitted by incoming particles 
will be converted into electron/positron pairs via the pair-creation process with the soft 
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X-rays and/or magnetic field. Those pairs are mainly produced at the convex side of the 
magnetic field lines, on which the primary 7-rays were emitted. New born pairs will loose 
their energy via the synchrotron radiation and/or the resonant Compton scattering, in which 
the Compton scattering is occurred in the resonant energy T^^Ex ~ tieBd{r)/mf.c with Fgi 
being Lorentz factor of the pairs. Those photons could also make new pairs. This cascade 
process will produce lots of low-energy electron/positron pairs near the stellar surface. The 
incoming pairs will carry the energy of Lg± ~ XLj. ~ 4 x lO'^^erg s~^. 

The incoming pairs eventually reach the stellar surface and heat up the stellar surface. 

The temperature of the heated platelet may be kB{Le±/i'^asBY^^^ 2(Le±/10^^erg s~^Y^^ 
(f/lO^cm)^!/^ i^gY^ ^j^j^ 

being the Boltzmann constant and asB Stephan-Boltzmann 
constant. These hard X-rays emitted from the heated platelet may be too high to escape 
from the resonant Compton scattering near the stellar surface. The scattered X-rays will be 
redistributed on the stellar surface, and emerge as the soft X-ray emission with a temperature 
of kTs ~ k{Le±/AnRlasBy^^ ~ 0.6{Le±/W^^eTg s'^)An~^/^ keV, were Afi is the solid 
angle of the emission. It is also possible that because some fraction of the incoming pairs 
will be created outside the oscillating magnetic field lines that connect to the cracked platelet, 
the typical size of the heated surface by the incoming particles is larger than i. In the present 
scenario, hence, the temperature of the intrinsic emissions from the heated surface by the 
incoming particles will be kT^ ~ 0.5keV. 

For "outwardly" propagating GeV 7-rays produced near but below r ~ 5 x 10^ cm, 
the cascade process will not develop to reduce the photon energy to ~MeV, and stop after 
producing the first or second generation of pairs. For example, the first generation of pairs 
produced by the magnetic pair-creation process will emit the synchrotron photons with an 
energy ~ ~ lOOMeV, where is the primary GeV 7-ray produced by the curvature 

radiation (Cheng et al. 1998). The outwardly propagating photons with ~ lOOMeV will 
escape from the magnetic pair-creation processes if they are produced beyond ~ 10^ cm 
(c.f. equation ([T])). Between rd ~ lO^cm and Tc ~ 5 x lO^cm, hence, outwardly propagating 
~ 100 MeV photons emitted by the first geration pairs will escape from the magnetosphere 
of the magnetars, and will be contribute to the observed spectrum. We expect that the 
luminosity of this component will be LiooA'/eV ~ AL^, where A ~ (rc — rd)/ {5 x lO^cm) ~ 0.1. 
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3.4 Pulse profiles 

In this section, we discuss the expected pulse profiles of the 7-ray emissions. We refer the 
previous studies (e.g. Takata, Chang & Cheng 2007) for the calculation of the pulse profile. 
We apply the vacuum rotating dipole field as the magnetic field in the magnetosphere. We 
describe the polar angle of the magnetic field lines at the stellar surface with a = 6{(j))/6p{(j)), 
where 6p is the polar angle of the last-open field line and (p is the magnetic azimuth. The 
open and closed magnetic field lines correspond to a < 1 and a > 1, respectively. Because the 
polar cap radius is only ~1% of the stellar radius, it is likely that the cracking platelet by the 
magnetic field is occurred at the closed field line region. As we discuss in section [21 the GeV 
7-rays can escape from the pair-creation process if they are produce beyond r ~ 5 x lO^cm 
from the stellar surface. In this section, hence, we apply a = 16, in which the magnetic field 
lines extend up to r ~ 5 x 10^ cm, as the main emission region. The width of the emitting 
magnetic flux tube in the azimuthal direction is assumed to be 50 ~ 1 radian, because the 



typical size of cracked platelet is RsOcScj) I lO^cm, where 9c ~ \J Rs/tc ~ 0.1 radian is 

the polar angle of the magnetic field lines that extend beyond r ~ Tc ~ 5 x 10 ''cm. 

We assume a constant emissivity along the magnetic field line, and at each point we 
express the opening angle of the 7-ray cone as 



which represents the effect of the oscillation of the magnetic field lines, on which the Alfven 
wave propagates (c.f. equation (|6])). The direction of the centre of the radiation cone coin- 
cides with the direction of the background vacuum dipole field. Because the Alfven wave 
propagates along the closed magnetic field lines, the wave may bounce many times between 
the footprints of the field lines. The crossing time scale of the Alfven wave, Tc ~ lO^cm/V^ ~ 
0.01s with Va ~ c being the speed of the Alfven wave, is much shorter than the rotation 
period. To calculate the pulse profile, hence, we take into account the 7-ray emissions from 
both particles migrating from north to south poles and from south to north poles. The 
inwardly propagating 7-rays may be absorbed by the pair-creation process with the back- 
ground X-rays and/or the magnetic field if they pass through near the stellar surface. To take 
into account this effect, we ignore the contribution from the inwardly propagating 7-rays 
that pass through the region r < 5 x lO'' cm. 

The predicted light curves for the inclination angle a = 60° are summarized in Figures [2] 
and [3], which are the sky-photon-mapping and the calculated light curves, respectively. We 





(17) 
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assume that the azimuthal angle of the centre of the emitting magnetic flux tube, which has 
a width (50 ~ 1 radian, is 0c = 0°,90°, 180° and 270°. In Figure [21 the whiteness refers the 
emissivity of the radiation on [C,, $)-plane, where C, and $ are the Earth viewing angle and 
the rotation phase, respectively. The brighter region (e.g. $ ~ 0.2 and ^ ~ 100° in the panel 
for 0c = 90°) corresponds to the rotation phase at which both radiations produced by the 
particles migrating from north to south poles and from south to north poles contribute to 
the observed emissions, while the dark region (e.g. $ ~ 0.75 and C, ~ 80° in the panel for 
0c = 90°), both emissions passes through the region r ^ 5 x 10^ cm and do not contribute 
to the observed emissions. 

In Figure[3l from top to bottom, the calculated light curves are for C, = 30°, 60°, 120° and 
150°, respectively. We can see in Figure [3] that the calculated light curves have in general 
broad peak in one-rotation period with sometimes sharp peaks (e.g. 0c = 90° and ^ = 120°). 
For the emission from the closed field lines, the special relativistic effects can be ignored, 
while it becomes important to explain the sharp and narrow pulse profiles of canonical 7-ray 
pulsars (Romani & Yadigaroglu 1995), in which the observed emissions are produced on the 
open field line region. In the present scenario, on the other hand, the sharp peak appears if 
the radiations from both kinds of particles migrating toward south pole and toward north 
pole are observed simultaneously. 

3.5 Spectrum 

To calculate the typical spectrum, we assume that the acceleration region extends beyond 
r ^ 10'' cm and all of GeV 7-rays emitted above Tc = 5 x 10^ cm can escape from the pair- 
creation process. Between rd ~ lO^cm and Tc, furthermore, all synchrotron photons emitted 
by the outwardly migrating pairs, which were produced by outwardly propagating 7-rays, 
can escape from the pair-creation processes, as discussed in section 13.31 The spectrum of the 
curvature radiation at each calculation grid is expressed as 



where -ft's/s is the modified Bessel function of order 5/3. In addition, ja is the current density 
and A is the cross section of the grid perpendicular to the dipole field. In this section, we 




(18) 



Fix) = X 
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assume the Goldreich- Julian value, ja{f^) = ^Bii{r) /27t, for the current density. We assume 
that the Lorentz factor of the accelerated particles is given by equation ( fT6|) . 

The synchrotron radiation of the outwardly migrating pairs produced between Vd and 
Tc (c.f. section 13731) may be described by 

Ps{E,) = --^ / — £F C dEe, 19 

where ( = Eg/Esyn with Egyn = STl^ehB ii{r) sin Op/inmeC and 9p the pitch angle. The 
distribution of the pairs under the steady condition may be described as 



dEe ^ E, 



^e±,min ^ ^ e± ^ E^jlmeC?, (20) 

■'syn 

where Esyn = 2e'^B'^ sin^ 9pTlj_/3m'^c^ is the energy loss rate of the synchrotron radiation. 
The minimum Lorentz factor of the pairs is re±,mm ~ 1/ sin6'p ~ rrieC^ B,/ {xE-yBd{r)) . 

Figure H] summarizes typical spectra predicted by the present scenario. In Figure HI 
we have used a = 60° for the inclination angle, a = 16 for the magnetic surface of the 
typical emission region, and (pc = 0° for the azimuthal angle of the centre of oscillating 
magnetic flux tube. In Figure HI the right and left panels show the dependency of the 
spectra of the curvature radiation on the released energy (Etot) and on the global dipole field 
{Bii{Rs)), respectively. The results are for ^2 = 1 s^^ and i = lO^cm. We find in Figure H] 
that the luminosity is proportional to the released total energy because oc 5B^{Rs) oc 
i^tot; as indicated by equation f|TT]) . As we can see in Figure HI the typical energy of the 
curvature spectrum increases with the released energy, that is, Ec oc 6B'^/\Rs) oc Efl^^ (c.f. 
equation (fT4|) ). In the right panel of Figure HJ we see that the typical energy of the curvature 
radiation decreases with increase of the magnetic field. This is because the displacement of 
the footprints of the oscillating magnetic field lines becomes smaller for stronger background 
magnetic field, as equation f lTOj) indicates. The decrease of the displacement implies the 
decrease of the potential drop (equation ([9])) and hence the decrease of the typical energy 
of the curvature radiation, Ec oc B^^^^{Rs) (equation [T^ . 

In Figure O we present the predicted spectrum in wide energy band using the parameters 
of AXP IE 2259+586 {Q ~ O.ds'^ and Bd{Rs) ~ lO^^G). AXP IE 22459+586 has shown 
frequent glitches, which sometimes accompany the X-ray outbursts. Furthermore, the GeV 7- 
ray emissions in the direction of CTB 109, which will be associated with AXP IE 22459+586, 
were founded in the Fermi data (Castro et al. 2012). Although the origin from SNR has 
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been suggested, the possibility of the emissions from AXP IE 22459+586 has not been ruled 
out yet. 

The model spectrum in Figure [5] is result for Etot = 2 x lO^^erg, i = 5 x 10^ cm, (pc = 0° 
and a = 13. The solid line and dashed line are spectra of the curvature radiation above 
r ^ 5 X lO^cm and between lO^cm ^ r ^ 5 x lO^cm, respectively. The dotted line represents 
the synchrotron spectrum from the first generation of pairs, if all outgoing photons above 
IGeV emitted between lO^cm ^ r ^ 5 x lO^cm (dashed line) are absorbed by the magnetic 
field (c.f. section [3l3|) . The measured spectrum is represented by the filled circles of Figure |5l 
We can seen in Figure [5] that if the observed GeV emissions would be originated from the 
magnetosphere, the observed flux level could be explained by the present scenario with the 
total released energy of Etot = 2 x lO'^^erg and the typical size of the cracking platelet of 
£ ~ 5 X lO^cm, respectively. However, the present model can reproduce only emissions below 
~10 GeV, as Figure [5] shows. Above 10 GeV, the present model expects that the emissions 
from SNR dominate the magnetospheric emissions. 

The synchrotron radiation near ~ 5 x 10^ cm (dotted line) can extend down to hard X- 
ray bands, but its flux level is well below the persistent X-ray emissions measured by RXTE 
(filled triangles, Kuiper et al. 2006). Hence, it will be difficult to observe the synchrotron 
radiation in the outer magnetosphere. 

4 DISCUSSION 

We have discussed possible scenario for the pulsed GeV 7-ray radiation in outer magneto- 
sphere of the magnetar. The GeV 7-rays will escape from the pair-creation process if the 
emission process is occurred beyond Tc ~ 5 x 10^ cm. In the present scenario, the Alfven 
wave carries the magnetic energy released by the crust cracking of the magnetic field into 
outer magnetosphere r ^ along the background magnetic field. The oscillation of the 
magnetic field induces the available potential drop 5$p ~ 10^^ Volts, which can accelerate 
the electrons and/or positrons to the Lorentz factor F ~ 10''. The curvature radiation at the 
outer magnetosphere can produce GeV 7-rays. 

The pulsed GeV 7-ray radiation from the magnetars have not been reported yet, although 
the predicted luminosity ~ 10^^ erg s~^ may be large enough to detect pulsed GeV 7-rays 
by the Fermi telescope. However, several reasons can be raised to explain the non-detection 
of the pulsed GeV emissions from the magnetars. First, the magnetars are in general located 
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at the Galactic plane, the background radiation may prevent the detection of the pulsed 
radiation. Second, the magnetars have shown frequent glitches that are sudden changes in 
frequency and/or frequency derivative (Igdem et al. (2012)). Hence the timing parameters 
of the magnetars are very unstable, which makes even harder to detect the pulsed period 
in the Fermi data. Third, it has become clear that the magnetars are associated with 
the SNRs (Allen & Horvath, 2004; Gaensler et al. 2005; Halpern & Gotthelf 2010). It is 
possible that 7-ray emissions from SNRs dominate the pulsed emissions in the data. Recent 
Fermi observations have revealed properties of the GeV radiations from SNRs (Abdo et 
al. 2009, 2010, 2011). For the young SNRs, the accelerated electrons at the forward shock 
may produce the GeV radiations through the inverse- Compton process, while for middle 
age SNRs, the interaction between the accelerated protons and ambient molecular clouds 
produces the bright GeV radiations via the ir^ decay process. 

We note that the typical luminosity of GeV radiation of SNRs is order of 10^^~^^ erg s~^, 
which is same order of magnitude predicted by the present magnetosphereic emission model. 
Hence, it would be possible that the GeV emissions in the direction of magnetars are com- 
posed of the emissions from SNRs and magnetospheres. The pulsed radiation predicted by 
the present scenario will change its luminosity level ~ 10'^^~^^erg s~^ at a temporal scale 
of years after the energy injection into the magnetosphere, while the SNR's emission will be 
stable. Hence, a temporal behavior of the observed GeV emissions will discriminate between 
the two components. 

Finally we note that it has been proposed that the glitch of normal pulsars is caused 
by unpinning of the super- fluid vortices from the lattice (Alpar 2001). The energy released 
by the glitch is deposited at inner crust, and most of energy is dissipated to heat up the 
entire star (Tang & Cheng 2001). The increases of the surface X-ray emission may affect 
the gamma-ray emission because of the photon-photon pair-creation process. With a typical 
released energy ~ lO^^erg, the glitch can increase surface X-ray emissions by a factor of 2-3 
if the core temperature is ~ lO^K for old pulsars or only a few percent if ~ lO'^K for 
young pulsars (Tang & Cheng 2001). Hence, it is expected such small change of the surface 
emission does not change much the 7-ray fluxes for young pulsars, e.g. Crab and Vela. For 
the magnetar case, on the other hand, the cracking of the outer crust could trigger the 
glitches that accompany the radiative outburst. A part of the released energy is carried into 
the magnetosphere via Alfven waves, and may be used to accelerate the particles, which 
results in the GeV 7-ray emissions. 
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Figure 1. Schematic view for the pulsed GeV 7-ray emissions from magnetars. The figure represents snap shot at which the 
Alfvcn wave propagates from the north pole to south pole. The GeV 7-rays produced above rc ~ 5 X 10^ cm can escape from 
the pair-creation process. The created pairs below eventually heat up the stellar surface and will contribute to the soft X-ray 
emissions from the stellar surface. Some of outwardly propagating 100 MeV photons emitted via the synchrotron radiation of 
the first generation of pairs will also escape from the pair-creation process, if they are produced beyond rd ~ lO'^cm. 
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APPENDIX A: PROPAGATION OF ALFVEN WAVE 

Within the frame work of a simple geometry of the co- rotating magnetosphere (c.f. Figured]), 
we discuss the Alfven wave propagating along the magnetic field line. We define the z-axis 
as the spin axis, and assume an uniform magnetic field with x and z components. 

Bo = B,e, + B,e,. (Al) 

The co-rotation velocity is expressed as 
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Figure 5. The spectrum of AXP IE 2259+586. The soUd and dashed hnes are curvature radiation above and below r = 
5 X 10^ cm, respectively. The results are for Etot = 2 X 10*^erg and ^ = 5 X lO^cm. The dotted line is spectrum of the 
synchrotron radiation from the first generation of the pairs produced lO'^cm ^ r ^ 5 X lO^cm, if all curvature photons above 
IGeV (dashed line) are absorbed by magnetic field. The data arc taken from Castro et al. (2012) and from Wu et al. (2012) 
for the Fermi and from Kuipcr ct al. (2006) for the RXTE, respectively. 
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Figure 1. Geometry of the Alfven wave propagating along the magnetic field line. The unperturbed fluid rotates around the 
spin axis, v = v^e^ + VyCy. 



VQ = V^e^ + VyBy 



(A2) 



where = —^y and Vy = Qx, respectively. We will assume a uniform perturbation of the 
magnetic field 
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5B = 6Bey, (A3) 

where SB is constant in both x and time t. The hnear differential equations for the equation 
of motion and the Maxwell equations are written as 

6F, (A4) 



P 



dt 

V ■5E = Andpe, (A5) 
V-6B = 0, (A6) 
— (A7, 
and 

where the perturbation of the Lorentz forces can be written down as 

5F = SpeE + pJE + ^{5j xBo + jx 5B), (A9) 

where pe and j = p^Vo are background charge density and current, respectively. The back- 
ground charge density, pe, is described by the Goldreich- Julian charge density, that is, 
Pe = V ■ E/Air = —V ■ {vq X Bo)/4vr ~ —VLBz/2ttc. The remaining equation is given 
by the frozen in condition 

5E = --{5v X Bq + vqX 5B). (AlO) 

c 

We may seek a solution in the form of the plane wave propagating along the background 
magnetic field; 

f (xey.\;){i[k{hxX + hzz) — ut]}, (All) 

where = Bj./ B, = B^/B, uj is the wave frequency and k is the magnitude of the wave 
number. As this take place in the set of the linear equations (lA4p - (lA9p . we can obtain the 
electric field as 

5E^ = -b,{-P^ + b,l3,)5B, 6Ey = 0, 6E, = -%I3^ + /3.6^)55, (A12) 
and the current as 

53x = -t-ribz5B - I3^5E^), 6jy = 0, = i—{bJB + I3^5E,), (A13) 
47r 47r 

where = u / kc is the phase velocity in units of the speed of light and (3^ = v^/c. 
The dispersion relation is given by 

(fe./3. - - ^ + b'j'^A - - = 0, (A14) 
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where /Sa = va/c with va = Bo/y/Airp being the Alfven velocity. Usually, the second term 
in the left hand side can be negligible because Q/kc ~ Q/u ^ 1. In the limit of f3x Pa 
and f3y = Vy/c ^ /3a, the phase velocity becomes 

(A15) 

which is the traditional dispersion relation for the Alfven wave, and we obtain = 1 for 
Pa » 1 and = (3a for (Sa < 1. 

The electric field parallel {SE\\) and perpendicular {6E±) to the background magnetic 
field can be described by 

6Eii = \6E ■ Bo\ = \h(3JB\ (A16) 
and 

6E^ = \6E X Bo\ = 1/3/51, (A17) 
respectively. For the current, we obtain 

= -i^KPJB, (A18) 

and 

6j^ = t^{l-Pl)5B, (A19) 

respectively. The important result is that the electric field and current along the magnetic 
field are exited due to the rotation of the unperturbed matter, and hence the electrons 
and/or positrons can be accelerated along the magnetic field line. 
In the limit of /3a ^ 1 (/3(<> = 1), we find that 



(<5£;|,, 5E^) ~ -6B{hA. 1). (A20) 

Hence the magnitude of the induced electric field is in order of magnitude of the perturbed 
magnetic field. The current becomes 

~ i^5B{^hJ^, f3f), (A21) 

implying the parallel current dominates the perpendicular current. 

In the limit of P^^y <^ /3a ^ 1 {P^ = /3a), the electric field and the current become 

6E^)r^-6B{bA,P^) (A22) 

and 

~ z£55(-6./3., /3a^), (A23) 
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respectively. The magnitude of induced electric field is much smaller than that of the per- 
turbed magnetic field. 

Finally, if (3 a ~ (3x,y <^ 1, we will see the phase velocity 

/3^~Mx±/3a. (A24) 



